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Abstract 



We study the local time distribution of a Brownian particle diffusing along the links on a 
graph. In particular, we derive an analytic expression of its Laplace transform in terms of the 
Green's function on the graph. We show that the asymptotic behavior of this distribution 
has non-Gaussian tails characterized by a nontrivial large deviation function. 



Graphs are ubiquitous and fascinating objects Q. In equilibrium statistical physics, the 
study of various model systems on graphs have provided deeper insights into how thermodynamic 
properties depend on the geometry of the underlying graph or network. For example, the study of 
a spin system on a complete graph provides an understanding of the thermodynamical properties 
at a mean field level. These studies on graphs are often nontrivial, the famous example being 
the Sherrington-Kirkpatrick model of Ising spin glasses defined on a complete graph [^. More 
recently, models of nonequilibrium statistical physics such as the abelian sandpile model have 
also been studied on graphs providing deep insights to the dynamics Q. The simplest dynamical 
example that has evoked enormous interest in both the physics and the mathematics community 
is the study of random walks on graphs, where a particle hops from one vertex to another 
provided they are connected by an edge. This has led to the study of the spectral properties 
of the discrete Laplacian operator on both regular Q and random graphs. All these systems 
mentioned above share one common property: the links between vertices on the graph do not play 
an active role in the actual physical process, their only purpose is to provide just a connection 
between two vertices. For example, the Euclidean length of a link is completely irrelevant and 
the physical properties only depend on parameters such as the coordination number of vertices 
etc. Due to this passive nature of the links, we refer to such networks as 'passive-link' networks. 

On the other hand, there exist 'active-link' networks where the links participate directly in 
the actual physical process. Examples include the networks for supplying household utilities 
such as electricity, water, telephone etc, the network of pipelines carrying oil and natural gases, 
the network of blood vessels in a living organism and many others. Quantum transport through 
mesoscopic networks provides yet another example of such 'active-link' network which has been 
studied experimentally and theoretically fs], Unlike their passive counterparts, the properties 
of physical observables on 'active-link' networks do depend on the Euclidean lengths of the links 



in addition to other parameters of the graph. Given the fact that the study of random walks 
on 'passive-hnk' networks has evoked so much interest and found numerous apphcations, it 
is natural to extend these studies on 'active-link' networks. In contrast to the 'passive-link' 
networks, the random walker does not hop between vertices of an 'active-link' network but 
actually undergoes continuous time Brownian motion along the links. A biological example of 
such a diffusion process on an 'active-link' network is the spread of an infectious virus or bacteria 
along the blood vessels in a living organism. 

In this Letter we study analytically the properties of a particular physical observable asso- 
ciated with the diffusion process on an 'active-link' graph, namely the probability distribution 
of the local time spent by the Brownian particle at a given point on the graph. This quantity 
characterizes the amount of time T spent by the particle between and t in the vicinity of a 
point. The corresponding quantity for a discrete random walk on a 'passive-link' network is the 
widely studied 'number of returns' to a given vertex. In the biological context mentioned above, 
the local time T denotes the time spent by the diffusing virus in the vicinity of a point (e.g, 
near the brain or the lungs) within its own lifetime t and hence is a measure of the damage that 
the virus can cause at a particular place in the network of blood vessels. Thus the study of the 
statistical properties of the local time in this context is important for medical purposes. 

For a fixed total time t, the local time T at any given point on the graph is clearly a random 
variable taking a different value for each history of the diffusion process. In this Letter we study 
analytically the probability distribution P{T,t) of the local time and show that for any generic 
graph where the total lengths of the links is finite, the distribution has the generic asymptotic 
behavior, P(T,t) ~ exp [—t^{T/t)] in the scaling limit T — > oo, t — > oo but keeping the ratio 
T/t fixed. The function <I>(x) is a large deviation function that characterize the non-Gaussian 
tails of this distribution. We provide a general formula for this large deviation function <I>(x) 
and calculate it explicitly in few specific examples. 

Let us begin by establishing our notation. 

Consider a graph Q consisting of a set of V vertices, labelled from to V — 1 and linked by B 
links of finite lengths. The coordination of vertex a is denoted by rria, therefore J2a=o ~ 

Each link [af3] of length lap is identified with an interval [0, lais] of R. We denote by Xaf^ the 
coordinate on the link [a(3] starting from vertex a . Unless otherwise stated the total length of 
the graph / = J2[ai3] is assumed to be finite. 

Now, let us consider a Brownian particle starting at t = from some point O. Without 
loss of generality we may assume that this point is a vertex (a = 0) of the graph. As time 
evolves the particle will diffuse along the links and explore the whole graph. If the diffusion is 
recurrent (which is, in particular, the case for a finite graph), it will revisit infinitely often the 
initial vertex. The aim of this letter is to compute the probability distribution of the total time 
T spent at O between and t. If x{t) denotes the position of the particle at time t, the local 
time T spent by this particle in an infinitesimal neighbourhood of O may be defined as 

f5{x{t')) dt', (1) 
Jo 

where x = is the location of O. 

Let -P(T, t) be the probability distribution of T. This distribution was studied recently in 
the one dimensional Sinai type models and it was shown that the characteristic function of the 
local time distribution can be expressed in terms of the evolution of a quantum Hamiltonian 
by using a path integral or equivalently a backward Fokker-Planck approach |^]. This method 
can be easily extended to the present case where one dimensional lines connect up to form the 
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'active-link' network. We find 
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/ dT P{T,t) e-P^ = / dx 
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where E denotes the expectation with respect to all Brownian paths over time t each starting 
from O and H is a quantum Hamiltonian defined below. A further Laplace transform in t gives 
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where G{x) is the resolvent. On each link [a(3], the Hamiltonian H acts as the one dimensional 
Laplace operator — (= )• Moreover, the behavior of the resolvent G{x) has to be 

a/3 

specified in the neighbourhood of all the vertices. 



Consider indeed some vertex a and its nearest neighbours i = 1, 2, . . . , rua- Suppose that the 
Brownian particle reaches a. It will leave this vertex along the link [aPi] with some probability 
that we may choose arbitrarily. For the sake of simplicity, we make the homogeneity 
assumption pa/s^ = l/m^. With this assumption, the resolvent G{x) is shown to be continuous 
at all the vertices. We will denote by Gq, its value on vertex a. 

If a 7^ 0, probability conservation implies 



E 

1=1 



dG 



dx. 
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In addition, in an infinitesimal neighbourhood of the starting vertex O, G must satisfy 

1 



-A + -f + p6]G = 6. 



(4) 



(5) 



This equation is a generalization of the one dimensional case studied in Ref. |Q . 

Denoting by one of the nearest neighbours of O, spatial integration in this neighbourhood 
gives 

1 ^ dG 

- E A + pGo = I. 6 

Let us now show that all the derivatives of G appearing in the above equations can be expressed 
in terms of the Gq,'s. On the link [a/3], G{xai3) satisfies 



1 d^ \ 

2 ^ ^ ) GiXafu) = 0, 



2 dx 



whose solution is 



It follows that 



^[Xal3) — ■ , /T^, 1- 
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Cap 
Sap 



27 COth \/2jlap = C/3o 
27 

Spa- 



sinh ^/2rflap 
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The relation (P) allows us to write Eqs. in a matrix form 



M{p) G = L, (12) 
where M(p) is a symmetric {V x y) matrix with elements 



mo 



Moo = ^co;,, +2p (13) 



m„ 



M„„ = ^c„ft if a/0 (14) 



=1 



Mx^ = -sxf, if [Xfj] is a bond (15) 
= otherwise. (16) 

G and L are two (y x 1) matrices whose elements are, respectively, Ga and Lq, = 25q,o- Inverting 
Eq. (0) we get 

G„ = 2 fM(p)-i) ^ . (17) 



' qO 

Substituting the results from Eqs. (|^) and (|l^ in the original expression in Eq. (|3|) and carrying 
out a few elementary manipulations on determinants, we get a closed form expression for the 
(double) Laplace transform, 

1 detM(O) 

7detM(p)- ^ ^ 

To proceed further, let us denote by Ga the value of Ga for p = 0. Due to the special form of 
the matrix M{p), we may write 

det M(p) = (1 + p Go) det M(0) (19) 

and recast ( [l8|) in the form 

roo roc 1 1 

/ dTe-P'^ dt e-^' PiT,t) = - (20) 

^0 Jo 7 1 + p Go 

Since Go is independent of p, we may invert the Laplace transform with respect to p and get 
the relationship 

r dt e-^' P{T, t) = ^e-^^W (21) 
Jo 7 

with 

A(7) = i. (22) 
Go 



Eq. ( [2l| ) is the central result of this Letter. It provides a generalization to graphs of the 
Levy formula for one dimensional diffusion processes ||l^. It is worthwhile to notice that Eqs. 
(1^.21, 22|) are valid for all kind of graphs irrespective of whether they are finite or not. The 
simplicity of this result calls for an alternative derivation based on paths integrals. The starting 
point is a representation of the resolvant as a two point function of a scalar field theory where 
the field (j){x) is a scalar field defined on the links of the graph and satisfying suitable continuity 
properties ((P). 

1 



= Z{p) (23) 
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where 



and the action is 



Z{p) 



/Pc/>(/>(x)(/>(0)exp-5((/>) 
/P(/.exp-5((/>) 



5(0) 



dx(^(V(/.)2 + 7(/>2 



(24) 



(25) 



One first expands the exponential term in p both in the numerator and denominator and com- 
putes all contractions using Wick's theorem. Then one can resum the series and obtain after 
some algebra 



X 



1 



(26) 



Substituting this result in Eq. 
(|2^ 



and integrating over x, we get back the same result as in Eq. 



The computation of the local time distribution is therefore reduced to the calculation of the 
Green's function on the graph. In some cases it is quite useful to express A (7) as a ratio of two 
spectral determinants of the form 

^<^' - \Wy 

where 5(7) = det(— A + 7) is the spectral determinant of the graph. The numerator is com- 
puted with Neumann boundary conditions at all vertices and the denominator is computed with 
Dirichlet boundary conditions at O and Neumann at all other vertices. Note that such ratios 
of spectral determinants which arise in the context of scattering theory on graphs ||8|, p!^, 13 



have recently appeared in the mathematical litterature in relation with Dirichlet forms or 
Sturm-Liouville problems [^. 

In order to extract the behavior of the distribution P(T, t) from the general formula in Eq. 



(21), we need to compute the function A (7) explicitly for a given graph and then invert the 
Laplace transform. In general, it is not easy to invert this Laplace transform. However, in the 
asymptotic limit when both t and T large with their ratio fixed, it is possible to make progress. 
In this scaling limit, one expects, from generic considerations, the following scaling behavior 
§0, 

P(r,t) ~e-**(^/*), (28) 

where $(y) is a large deviation function. Substituting this expected scaling form on the left 
hand side of Eq. (|2^) and carrying out a steepest descent calculation valid for large t and T, 



one obtains the following Legendre transform. 



$(y) = max[-7 + y A(7)]. (29) 
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Thus the large deviation function can be computed from Eq. ( p9D provided one knows 
the function A(7) explicitly for a given graph. Examples where one can calculate explicitly 
the function A(7) are given below. Interestingly, however, some general features of the large 
deviation function $(y), such as its behaviour near the tails as well as near its minimum can be 
worked out for a generic graph even without the detailed knowledge of the function A (7), thus 
displaying certain universal features. 

First, let us consider finite graphs. One can show that <I>(y) decreases monotonically with 
increasing y in the range < y < 1//, achieves a minimum at y = l/l and then increases 
monotonically for y > l/l (recall that / is the total length of the graph). We consider the 
following three regions where the behavior of $(y) is quite generic. 
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i) Let us first consider the hmit 7 ^ cxo. Using Eqs. (|l3Hl|) and (||), we get 



m-rnoJ'l, (30) 



where rriQ is the coordination of the vertex O. Using this expression in Eq. (29) and maximizing, 
we get 

2 2 

$(y) w !!!:^JL when y^oo. (31) 
8 

The hmit y — > cxd means that we select infinitely long Brownian trajectories (recall that i — > cxd) 
for which the particle stays during almost all the time close to O. Since the particle does not 
explore the whole graph, the local time distribution will mainly depend on the characteristics 
of the graph in the vicinity of O, i.e. on the coordination rriQ. 



ii) By considering the first negative singularity, 70, of A(7) in Eq. (|22|) we can show that, for 
7 « 7o, A(7) ~ 0/(70 — 7) with a > 0. From Eq. ( [29| ) we deduce that 

«>(y) « |7o| when y ^ 0+ (32) 

where 70 is the lowest eigenvalue of the Laplacian on the graph with Dirichlet boundary condi- 
tions in O. 



iii) Finally, it turns out that 7 — > limit corresponds, via Eq. (|29|), to y — > In this limit 



I f iV , 1 



(34) 



The positive constant C is obtained by a small 7 expansion in Eq. (jl^) with the result: 

detg 



detN 



where the matrix is obtained by removing the 1st line and 1st column in the matrix M(0) 
and taking the limit 7^0. The matrix Q is obtained by taking the limit 7 ^ in the matrix 
M(0) and replacing the 1st line and 1st column by the elements Qoo = ~ (1/3) [a/3] ^a/? 

Qoa = QaO = ET=1 lap, ^^V O / 0. 



Note that the value l/l corresponds to the mean local time {T)/t and the constant C is 
proportional to the variance of the local time. This result, which also follows from a cen- 
tral limit theorem, indicates a Gaussian distribution for the local time near its mean value, 
P{T,t) ~ exp [— (T — t//)^/2(T^], where the variance, o"^ = Ct/2 for large t. This immediately 



gives a large deviation function ^{y) as in Eq.(33). Note, however, that this central limit theO' 



rem holds only near y k, l/l but fails in regions away from y = l/l. In particular, near the tails 
the distribution becomes non-Gaussian as refiected in the form of $(y) near y — > and y ^ oo 
above. For one dimensional diffusion processes, these characteristics of large deviation functions 
were identified not just for the local time but for other functionals of Brownian motion as well 
i0. 

As an example, consider a complete graph Kn with n vertices where there is a link of length d 
between any pair of vertices (see fig.l). In this case, putting p = in Eq. (p^, evaluating the 
inverse of the M matrix explicitly and using this result in Eq. (|22|), one obtains (Q) 



w ^ 1 [cosh(dV2^) + cos((^)] 

2 cos((/>) [(cosh((i^) + cos(,/.) - 1)] '''' 



2^ ' 



(35) 
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Figure 1: The complete graph K^. 



where cos{<f>) = In particular, a small 7 expansion in (|2^ ) leads for the constant C to the 
value C = 16(n + 1/n — ll/6)/(n(n — 1))^, in agreement with Eq. p^). 



Another interesting example, which was recently used in the context of quantum chaos J16[ , 
1 17 1 is the case of a star graph consisting of n links of length d attached to the vertex O (fig 2). 




n-1 



Figure 2: The hydra with n legs. 

In this case one obtains 

A(7) = ny|tanh(dv^) (36) 

One can also compute explicitly the tails of the large deviation function and thus check that 
they are in agreement with the general results Eqs. ( |3l[]3^ ). Clearly, in the latter cases ii) and 
iii), the whole graph is explored and the distribution depends, in a rather nontrivial way, on the 
details of the graph. 

The previous formulae can be easily generalised to infinite graphs. The example of a star 
graph with infinite legs is interesting since it was already used in a probabilistic setting to derive 
multidimensional extensions of the arc sine law ( ||T^ ) . Using Eq. (^) one can derive the exact 
probability distribution 



Tl 

P(T,t) = -=exp 



8t 



(37) 



The fact that the distribution is purely Gaussian depends obviously heavily on the symmetries 
of the graph. However the behaviour for small and large T is universal. For an arbitrary graph 
with n infinite legs attached in any vertex one can show that $(?/) is a monotonically increasing 
function for y > and also that 



i) 



2 2 

~ !!}oy_ ^^^^ y^oo (38) 
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(same result as before and for the same reasons) 
ii) 

^{y) - when y^O (39) 

8 

(the particle tries to explore the whole graph; so, it essentially "sees" the infinite legs). 
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